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Abstract 

We propose a model for neutrino mass generation based on both the tree-level seesaw 
mechanism with a single right-handed neutrino and one-loop radiative effects in a non- 
supersymmetric framework. The generated mass matrix is composed of two parts which 
have the same texture and produce neutrino mass eigenvalues and mixing suitable for 
the explanation of neutrino oscillations. The model has a good CDM candidate which 
contributes to the radiative neutrino mass generation. The stability of the CDM candidate 
is ensured by Z 2 which is the residual symmetry of a spontaneously broken U(l)'. We 
discuss the values of U e ^ and also estimate the masses of the relevant fields to realize an 
appropriate abundance of the CDM. 



* e-mail : j ik@hep . s . kanazawa-u .ac.jp 
e-mail: suemat su@ hep . s . kanazawa- u . ac . j p 



1 



1 Introduction 



Recent experimental and observational results on neutrino masses [1] and cold dark mat- 
ter (CDM) [2] suggest that the standard model (SM) should be extended by introducing 
some neutral fields. A well studied candidate for the extension is the minimal super- 
symmetric SM (MSSM). Although the MSSM contains a good CDM candidate as the 
lightest superparticle (LSP) as long as the -R-parity is conserved, the parameter regions 
preferable for the explanation of the WMAP data are found to be strictly restricted in 
certain types of the MSSM [3]. Confronting these situations, it seems to be interesting 
to consider models in which we can explain these new features from the same origin in a 
non-supersymmetric extension of the SM: A certain symmetry related to the smallness of 
neutrino masses can guarantee the stability of a CDM candidate. Backgrounds that forth 
coming collider experiments like LHC may find signatures of such extended models make 
this kind of trials worthy enough at present stage. Several recent works have been done 
along this line [4]. 

In this paper we follow this line to propose an extension of a previously considered 
model by introducing a local U(l)' symmetry at TeV regions. As in the radiative mass 
generation models [5], we introduce an additional SU(2) doublet i] T = (i] + ,i] ) to the 
ordinary Higgs doublet H T = (H + ,H°). We also introduce a singlet whose vacuum 
expectation value breaks U(l)' symmetry spontaneously down to Z 2 which is responsible 
for the stability of the CDM candidate. This extension seems to remedy defects in the 
previous models that certain fine tunings are required for both the generation of small 
neutrino masses and the reconciliation between the CDM abundance and the constraints 
from lepton flavor violating processes. Based on such a model we calculate the value of 
the element U e3 of the MNS matrix and masses of the relevant fields which produce an 
appropriate abundance of the CDM. 

2 A model 

We consider a model with a similar symmetry to the model in [4, 6]. We extend it by 
introducing a singlet Higgs scalar <fi. This extension makes the model able to contain an 
additional U(l)' symmetry. In this paper we assume that this symmetry is leptophobic and 
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then leptons do not have its charge for simplicity. 1 The U(l)' charge for the ingredients 
of the model is shown in Table 1, in which fermions are assumed to be left-handed. 
Note that we need only two right-handed neutrinos Ni and N 2 to generate appropriate 
neutrino masses and mixings in a minimal case. Then the invariant Lagrangian relevant 
to the neutrino masses can be expressed as 

C m = E {h al L a HN 1 + h oa L ^N^ + ^M^ + ]-X<l>Nl + h.c., (1) 

where we assume that Yukawa couplings for charged leptons are diagonal. The most 
general invariant scalar potential up to dimension five may also be written as 

(f)(rfH) 2 + h.c." 

+ (m 2 H + \ 7 ^<f>)H^H + (m 2 + \ 8 ^<f>Wv + (2) 

We add a nonrenormalizable A6 term and a bare mass term for N\. The scalar potential 
(2) without the A6 term has an accidental U(l) symmetry, which forbids the one-loop 
contribution of the rj exchange diagram to neutrino masses. This symmetry is explicitly 
broken by the Yukawa interactions (1), so that terms like the A 6 term, i.e. (<ft '4>) n 4>(rf H) 2 , 
can be generated in high orders in perturbation theory in general. All of them contribute 
to radiative neutrino masses. 2 Here we do not ask the origin of the A 6 term. They might 
be supposed to be effective terms generated through some dynamics at an intermediate 
scale M*. We can check that there are no other dimension five operators invariant under 
the above mentioned symmetry in the scalar potential. 

As the model discussed in [4] , H plays the role of the ordinary doublet Higgs scalar in 
the SM but rj is assumed to obtain no vacuum expectation value (VEV). A singlet scalar 
is assumed to obtain a VEV, which breaks U(l)' down to Z 2 ( See Table 1). This VEV 
also gives the mass for N 2 through M N2 = A(0) and also yields an effective coupling for 
the A 6 term as \ 6 ((j))/M*. It can be small enough as long as ((f)) <C M* is satisfied. Thus, 
the masses of the real and imaginary parts of 77° are found to be almost degenerate. They 

1 We need to introduce some additional fermions to cancel gauge anomaly. Since such extensions will 

be done without changing the following results, we do not go further into this problem here. 

2 It turns out that one-loop corrections generating the \ 6 term, i.e. {(jy 4>)4>{rf H) 2 , vanish if the 

condition (5) discussed later is satisfied. 



+ X 4 (H^H)^rj)+X 5 (H^H) + -^- 
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Q a U a D a L a E a JVi N 2 H rj 



U(l)' 2q -2q -2q q -g -2g 



Z 2 +1 +1 +1 +1 +1 +1 -1 +1 -1 +1 



Table 1. Field contents and their charges. Z 2 is the residual symmetry of U(l)'. 



are expressed as M 2 



m2 + (A 4 + A 5 )( J f7°) 2 + A 8 (0) 



2 . In the model discussed in [4], the 



coupling constant of the term corresponding to this A 6 term is required to be extremely 
small to generate appropriate neutrino masses. This point is automatically improved by 
introducing the new U(l)' symmetry. 

3 Masses and mixings of neutrinos 

We find that there are two origins for the neutrino masses under these settings for the 
model. One is the ordinary seesaw mass induced by a right-handed neutrino Ni [7] and 
another is one-loop radiative mass mediated by the exchange of if and N 2 [5, 6]. These 
effects generate a mass matrix for three light neutrinos. It is expressed by 



Although two terms of M v may be characterized by different mass scales, the texture of 
both terms is the same as found in eq. (4). This type of the texture for neutrino mass 
matrix has been studied in [7, 8]. We neglect CP phases in the following discussion. 

Now we study eigenvalues and the mixing matrix for the neutrino mass matrix (3). We 
consider to diagonalize M v by using an orthogonal matrix U in such a way as U T M U U = 
diag(m 1; m 2 , m 3 ). If Yukawa couplings satisfy a condition 




(3) 



where v = (H°) and /i^ is defined by 




(a = 1,2). 



(4) 



h el h e2 + h^h^ + h Tl h T2 oc [fi^ , /i (2) ] = 0, 



(5) 
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jjS 1 ^ and fi^ can be simultaneously diagonalized. Since U can be analytically found in 
such cases, we confine ourselves to these interesting ones. We define a matrix U as 

/ 1 n n \ / pnc Q„ Gin 





cos 9 2 sin # 2 
^ — sin 9 2 cos # 2 J 



cos 63 sin 63 
1 
^ — sin 9\ cos 9 5 J 



(6) 



The first term of M v can be diagonalized by this U if the following conditions is satisfied: 



tan ^2 = -r— , tan# 3 = , 



(7) 



Then the mass eigenvalues for the first term of M u are obtained by using the following 
eigenvalues of //W: 

/i« g = diag(0, 0, h^ + h^ + hlJ. (8) 

We consider diagonalization of ^ next. At first, it should be noted that ^ is 
transformed by the same U. However, if the condition (5), which can be written as 



h e2 sin # 3 + {hya sin 9 2 + h r2 cos 9 2 ) cos 9 3 = 



(9) 



is satisfied, can be diagonalized by applying 
plemented by an additional one given by 



an 



orthogonal transformation UU3 sup- 



U 3 = 



cos9 1 sin9i 
— sin 9i cos 9\ 
1 



V 







(10) 



This additional transformation by U 3 does not affect the diagonalization of Conse- 
quently, both terms of M v can be simultaneously diagonalized by setting 

tan 9 2 tan 9 2 + 1 

tant^x = , (11) 

(tan 9 2 — tan 9 2 ) sin 9 3 

where we define 9 2 as tan# 2 = h^ 2 /h r2 . Finally, we obtain nonzero mass eigenvalues of 
the light neutrinos as 



tan 2 9 -\- 1 ~ A 

m 2 = AB - — ytt -(tan# 2 - tan# 2 ) 2 , m 3 = -(tan 2 9 2 + l)(tan 2 3 + 1), (12) 

tan u 2 + 1 2 

where A = 2h 2 rl v 2 /M* and B = (A 6 / 16tt 2 A) (h r2 /h Tl f I (M 2 2 /M 2 ). 
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Fig. 1 U e 3 as a function of B. Allowed regions are shown as the regions surrounded by red solid lines. 
Horizontal dashed lines stand for the present experimental upper bounds for | U e % | . 

Here we fix tan 62 = 1 which is supported by the data of the atmospheric neutrino and 
K2K experiment. CHOOZ experiments give the constraint on 83 such as | sin #3! < 0.22 
[9] . f f we use these conditions, the mixing matrix U = UU3 can be approximately written 

as 



U 



( 


COS^i 


sin Q\ 


sin 8^ 
V2 


\ 




sin 8\ 


cos 8\ 


1 






V2 


V2 


s/2 




V 


sin 61 


cos 81 


1 


/ 


V2 


V2 


v/2 



(13) 



Only two mass eigenvalues M v are nonzero and then we impose that squared mass 



differences required by the neutrino oscillation data satisfy = J Aml tm and m2 



Ar7ig 0l . Although there is a possibility that two nonzero eigenvalues have almost de- 



generate values such as y ^ m ltm an d their squared difference is given by Am^, we do 
not consider it since and have independent origins. We suppose Q\ = 9 so \, where 
^soi is a mixing angle relevant to the solar neutrino. Then we can determine 83 through 



eq. (11) by using tan0 so i, yAm^m, yAm s 2 ol and B. If we use neutrino oscillation data 
for these, we can find allowed regions of 63 as a function of B. This is shown in Fig.l, 
where we have used values of the measured neutrino oscillation parameters [10] 

Am s 2 ol = 8.0±g;| x 10" 5 eV 2 , Am 2 tm = (1.9 - 3.6) x 10" 3 eV 2 , 

tan 2 # sol = 0.45±°$. (14) 

This figure shows that B is restricted in narrow regions such as 0.03 < B < 0.1. 



6 



As an example, let us assume M v o/M N2 = 0.3-0.7 and then I(M%JM%>) = 0.1 - 1.3. 
In such cases h r2 /h Tl ~ 10(A/A 6 ) 1 / 2 should be satisfied. If we obtain more constraints on 
the relevant coupling constants, we may restrict the value of U e3 much more. Although 
U e3 takes a nonzero value for 0.03 < B < 0.05 and 0.08 < B < 0.1, U e3 = is also allowed 
for 0.03 < B < 0.08. The condition for the coupling constants can be easily satisfied 
even if we assume that coupling constants are 0(1). Therefore, the model needs no fine 
tuning to be consistent with all the present experimental data for neutrino oscillations. 
The effective mass m ee for the neutrinoless double beta decay takes the values in the range 
\m ee \ < 6.3 x 10" 3 eV. 

4 Relic abundance of a CDM candidate 

The lightest field with an odd Z 2 charge can be stable since an even charge is assigned 
to each SM content. If both the mass and the annihilation cross section of such a field 
have appropriate values, it can be a good CDM candidate as long as it is neutral. As 
found from Table 1, such candidates are N 2 and rf . Since they have a new U(l)' gauge 
interaction, their annihilation to quarks is considered to be dominantly mediated by this 
interaction. 3 If their annihilation is mediated only by the exchange of 77° or N 2 through 
Yukawa couplings as in the model discussed in [4], we cannot simultaneously explain, 
without fine tuning of coupling constants, both the observed value of the CDM abundance 
and the constraints coming from lepton flavor violating processes such as /i — * cy. Since 
U(l)' is supposed to be a generation independent gauge symmetry, we can easily escape 
this problem by assuming that the Yukawa couplings h a2 are small enough or both rf 
and N 2 are heavy enough. In the following study we consider the case that N 2 is lighter 
than rf. As seen in the last part of the previous section, this case is consistent with the 
present experimental bounds for U e3 without fine tuning. 

Now we estimate the relic abundance of N 2 and compare it with the CDM abundance 
obtained from the WMAP data. We suppose that possible annihilation processes N 2 N 2 — > 
// are dominantly mediated by the U(l)' gauge field. If it is expanded by relative velocity 

3 A role of U(l)' in annihilation of the CDM in supersymmetric models has been studied in [13]. 
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Fig. 2 Allowed regions by the WMAP data in the (Mz> , Mjv 2 ) plane. In both figures solid and dotted 
lines represent contours for fi^ft, 2 = 0.0945 and flN 2 h 2 = 0.1287. Green dotted lines stand for Mn 2 for 
typical values of A which are given in the text. In the right figure vertical dash-dotted lines represent 
lower bounds of M Z ' in case of \8\ = 10~ 2 , 5 x 10~ 3 , 10~ 3 . 

v between annihilating iVVs as av = a + bv 2 , the coefficients a and b are expressed as 

a = S c ^ Q ^ q (Trfe 5 ' b = ? + Q ^ )q TT^MW' (15) 



where (3 = y 1 — rnj/M^ and Cf=3 for quarks, s is the center of mass energy of collisions 
and q is the U(l)' charge of N 2 given in Table 1. The charge of the final state fermion / 
is defined as 

Q/V = QfR + QfLl Qf A = QfR ~ Qf L - ( 16 ) 

Using these quantities, the present relic abundance of N 2 can be estimated as [11], 



n i,2i M N2 n N2 
\L N2 n | 



Pcr/h? 



.76 x lfr 11 ^ 1/2 x F 



- (17) 



(a + 3b/x F ) GeV 

where g* enumerates the degrees of freedom of relativistic fields at the freeze-out temper- 
ature T F of iV 2 . T F is determined through the equation for a dimensionless parameter 
x F = M N JT F 

0.0955m pl M N2 (a + 6b/x F ) 
x F — In -— , (18) 

(g*x F ) 1 / 2 
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where m v \ is the Planck mass. If we fix the U(l)' charge of fields and its coupling constant 
g', we can estimate the present N 2 abundance using these formulas. Assuming a GUT 
relation g' = ^5/3^1 and q = 0.6 as an example, we calculate VLn 2 Ii 2 . The results are 
given in Fig. 2. 

In the left figure of Fig. 2 we plot favorable regions in the (M Z ,,M N2 ) plane, where 
VL^ 2 h 2 takes values in the range 0.0945 - 0.1285, which is required by the WMAP data. 
VL^ 2 h 2 has a valley in the parameter region of Fig. 2, and therefore the allowed regions 
appear as two narrow bands, each sandwiched by a solid line and a dashed line. Since 
Mn 2 and Mz> are induced by (0) and written as 

M N2 = \(<f>), M z , = 2V2g'q(<j ) }, (19) 

Mjv 2 is determined by Mz'- We plot this M^ 2 values by green dotted lines for A = 0.2 and 
0.7. The lower bounds of Mz> come from constraints for ZZ' mixing and direct search 
of Z'. H is assumed to have no U(l)' charge and then its VEV induces no ZZ 1 mixing. 
Moreover, since it is leptophobic, the constraints on Mz> obtained from its hadronic decay 
is rather weak. Thus, the lower bounds of M z > may be M z > ^ 450 GeV in the present 
model [14]. Taking account of this, Fig. 1 shows that this model can well explain the CDM 
abundance. Since A is included in the definition of B, values of 9% may be constrained by 
the mass of the CDM if we can obtain more informations on A 6 , h r2 /h Tl and M z >- 

Here we briefly discuss the relation to lepton flavor violating processes such as fx — > cy. 
As in the model of [15], fx — > is induced through the mediation of i]° and N 2 . Its 
branching ratio can be given by 

64vr(G F Mj)^|^ e2F2 (^)| ' 
F 2 (x) = 1 (1 - 6x + 3x 2 + 2x 3 - 6x 2 lux). (20) 

Taking account that 1/12 < F 2 (x) < 1/6 is satisfied in case of M N , 2 < M v o and imposing 
the present experimental upper bound B(fx — > erf < 1.2 x 10~ n , we find that M v o should 
satisfy 

M n o > (360 - 500) (^^j GeV. (21) 

Here we use the results of the previous section. Constraints coming from fx — > and the 
CDM abundance can be consistent for reasonable values of h r2 . Since N 2 annihilation due 



B(/x -»• erf 
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to an r]° exchange is ineffective for these values of couplings and masses [4], the results of 
the N 2 abundance given above is not affected by this process. 

Finally, it may be useful to refer to the cases of general U(l)'. In these cases a crucial 
condition for the mass of the U(l)' gauge field comes from the constraint for ZZ' mixing. 
A mass matrix for neutral gauge bosons can be expressed as 



\ 



(22) 



-g'jgl + gknv 2 2 g ' 2 q 2 (4(<j ) } 2 + v 2 ) j 

where qu and stand for the U(l)' charge of H and <fi. Since a ZZ 1 mixing angle 9 is 
known to be strongly suppressed [12], the magnitude of (0) should satisfy 

(0) ~ 2 (2q^\9W • (23) 
This condition gives a lower bound on both M z > and M N2 . In the right panel of Fig. 2 
we plot this bound in cases of \9\ = 1CT 2 ,5 x 1CT 3 , 10~ 3 , which are drawn by vertical 
dash-dotted lines. We also plot the values of M N2 for A = 0.2, 0.3, 0.6 and 0.9. They 
are drawn by green dotted lines. Although \9\ should be less than 10~ 3 , we may suppose 
larger values of \9\ in eq. (23) by extending the model without changing the results in the 
previous section. In fact, if the model has two Higgs doublets H u and which couple 
to up- and down-sectors respectively, off-diagonal elements of eq. (22) is proportional to 
g'iQHuiHu) 2 — lH d {H d ) 2 ) where qu ud expresses the U(l)' charge. Cancellation between 
these two contributions can make the ZZ' mixing smaller for the same value of (0). In 
such cases we can apply this effect by using larger \9\ values in eq. (23). In this figure 9 
values smaller than 10~ 3 should be understood based on this reasoning. 

On the other hand, the introduction of additional Higgs doublets may require us to 
take account of new final states for the N 2 annihilation induced by the Z' exchange. If 
N 2 is heavier than W ± , the final states should include gauge bosons and Higgs scalars 
such as W + W~, H®Hj, W ± H T , H + H~ and ZHf, where Hf is a mass eigenstate of the 
neutral Higgs. Since the annihilation to W + W~ is suppressed by the ZZ' mixing in the 
present model, important modes are expected to be H^H® and they may give the same 
order of contributions as the annihilation to // [11]. In order to take such effects into 
account without practicing tedious estimation of such processes, we show in the right 
figure of Fig. 2 an additional fl^h 2 contour which is obtained by using 5 x (crv)ff for 
cross section. It is drawn by blue lines. An original contour for the cross section (crv)fj 
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is drawn by red lines. 4 Since main parts of the cross section into these final states are 
expected to have the similar dependence on M z > and M N2 , this is considered to give good 
references for these cases. This figure also suggests that this kind of models can explain 
the CDM abundance even under the constraint for Z' physics. 

5 Summary 

We have studied neutrino masses and CDM abundance in a non-supersymmetric, but 
U(l)' symmetric model which is obtained from the SM by adding certain neutral fields. 
Neutrino masses are generated through both the seesaw mechanism with a single right- 
handed neutrino and the one-loop radiative effects. They induce the same texture which 
can realize favorable mass eigenvalues and mixing angles. One of the introduced neutral 
fields is stable due to an unbroken Z 2 symmetry which is the residual symmetry of the 
spontaneously broken U(l)'. Thus it can be a good CDM candidate. Since it has the U(l)' 
gauge interaction, the annihilation is dominantly mediated through this interaction. If 
this U(l)' symmetry is broken at a suitable scale, the present relic abundance of right- 
handed neutrinos can explain the WMAP result for the CDM abundance. This model 
suggests that two of the biggest questions in the SM, that is, neutrino masses and the 
CDM may be explained on the common basis of an extension of the SM. An interesting 
feature of the model is that the value of the third mixing angle 9 3 may be related to the 
mass of the CDM. The model may be examined through the search of the Z' and the 
additional Higgs doublet r\ at LHC. 
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4 In this calculation we use the U(l)' charge assignment for Higgs doublets, quarks and leptons such 
as Qh u = qH d = -2q, q Q = q L = 2q. 
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